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Abstract 

We analyse the redshift dependence of space number density of quasars assuming that 
they are the short-hved active stages of the massive galaxies and arise immediately after 
' the collapse of homogeneous central part of protogalaxy clouds. Obtained dependence 

O ■ fits the observational data ChaMP+CDF+ROSAT (Silverman et al. 2005) very well for 
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r~| , protogalaxy clouds of mass M ^ 8 - W'^ h-'MQ and ellipticity e < 0.4. The hfetime of 

Qh' bright X-ray AGNs or QSOs with Lx > 10^^'^ erg-s~^ in the range of energies 0.3 — 8 keV 

O ■ is TQso ^ 6-10^ years when the mass of supermassive black hole is Msmbh ~ 10^ Mq and 

the values of other quasar parameters are reasonable. The analysis and all calculations 
^ I were carried out in the framework of ACDM-model with parameters determined from 5- 

years WMAP, SNIa and large scale structure data (Komatsu et al. 2009). It is concluded, 
, that the halo model of galaxy formation in the ACDM cosmological model matches well 

^ I observational data on AGNs and QSOs number density coming from current optical and 

X-ray surveys. 
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^ Introduction 

The radio, optical and X-ray observational data indicate that co-moving space number density 
^ I of bright AGNs and QSOs is unmonotonous function of redshift: it increases up to 2; ~ 2.5, 
I declines and goes to zero after that |6l [HI [171 ESI EH [36]. Small angular diameters and huge 
luminosities of QSOs signify that their central power engines are supermassive black holes 
(SMBH) which efficiently transform the mass of accreting matter into radiation. About ten 
percents of mass can be converted into the high energy radiation by this mechanism instead 
of decimal parts of percent for the nucleosynthesis reactions providing the luminosity of stars. 
Observations using Hubble space telescope and the largest ground-based ones have revealed 
the host galaxies of QSOs in which they occupy central parts. High emissivity of quasars, their 
popularity in comparison with galaxies, genetic closeness to galaxy active nuclei support the 
idea that QSOs are short-lived active phases of galaxies with SMBH at their nuclei. So, the 
space number density of quasars can be described in the scenario of galaxy formation taking 
into account the peculiarities which make them QSOs. Really, in such approach the redshift 
distribution of QSOs space number density becomes qualitatively understood: the growth of 
QSOs number density in the redshift range z = (00, ~ 2.5] is caused by the rate of galactical 
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nuclei formation, the later decay in the redshift range 2; = [~ 2.5, 0] is due to the short-lived 
action of quasar engine (finite fuel resource in the vicinity of galaxy centrum and the large rate 
of its burning) as well as to the decreasing of galaxy birth-rate. In this paper we develop the 
semi-analytical approach proposed in the previous papers [21 EHl [IHl Ell [ID] in order to adjust 
quantitatively the prediction of concordance cosmological ACDM-model fT9] to observational 
data on bright QSOs number density given by [36j. 



1 Formation of bright QSOs in the halo model 

Dependence of the QSOs number density on redshift is determined by several factors: birth- 
rate of galaxies of relevant mass, probabilities of SMBH formation in their centra and the 
tempo of accretion of matter to them. The birth-rate of galaxies of given mass is determined 
completely by cosmological model and initial power spectrum of density perturbations. For 
the estimation of other factors the additional assumptions are needed. The high popularity 
of quasars suggests that the formation of SMBH is quite frequent phenomenon. The physical 
conditions in the centrum of forming galaxy are favourable for formation of the central SMBHs 
through the collapse of dark matter and hydrogen-helium gas [161 EH ISSl E2i El |22l E] as well as 
through the merging of remnants of supermassive stars [H [TSl [71 HI]. In both cases the time-lag 
between halo virialization and the birth of the bright quasar is expected to be short comparing 
to the cosmological time-scale, even at high redshifts. For appearance of the luminous quasar 
the one solar mass of gas per year must infall into SMBH. A few physical mechanisms of gas 
supplying to quasar engine are analysed in literature. For example, [28] suppose that accreation 
of gas into SMBH is caused by the tide collisions of close galaxies in the process of their merging. 
The semi-analytical models of the growth of SMBHs in the central parts of their host galaxies, 
developed by [261 123 EH IS], fit well the observational data on quasar luminosity function and 
two-point space correlation one for selected redshifts. 

In this paper we assume that bright QSOs appear in the process of galaxies formation 
and their luminous efficiency is provided by the collapsing gas from nearest vicinity of SMBH. 
The durations of both processes are essentially lower than cosmological time. So, quasars are 
supposed to be the short-lived active phases of the galaxies formation, which is well described by 
the halo model. The central density peak of protogalaxy cloud corresponding to the bottom of 
gravitational potential well is practically homogeneous and collapses first of all with formation 
of SMBH. The outer shells will collapse later and the relation between their co-moving radius 
R and the collapse time tcoi is given by simple formula following from the model of spherical 
collapse [ni ED]: _ 

M{R)-Mr 



■-D{1) = 6,{t 



R 

where M{R) = Att J p{r)r'^dr is the mass of the matter in the sphere with co-moving radius R, 


p(r) = p(l + S{r)), p is the mean matter density at moment tcoi, Mr = ^nR^'p, 6 is the initial 
amplitude of density perturbation, Sdtcoi) is the critical overdensity as function of collapse 
time, D{a) is growth factor to the present from high redshift relative to the growth factor from 
the same initial amplitude and initial redshift in the Einstein-de Sitter model, a is the scale 
factor, which equals 1 at current epoch. This formula is correct if mass of matter inside each 
shell is fixed. For close to centrum shells this condition is satisfied at the initial stage of SMBH 
formation because particles belonging to them have small target parameters. Outer shells with 
lower averaged density contrast will collapse later when hydrodynamical and dark matter halo 
virialization processes take place. Some part of matter will fall into a trap of SMBH providing 
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the quasar luminous efficiency. The quasar will be active some time after formation of galaxy 
and virialization of dark matter halo till the matter which serve the fuel for quasar engine will 
be depleted from the range of unstable orbits. The tide collision of merging galaxies can renew 
the quasar activity and this can be repeated many times during the galaxy life. We suppose 
that most luminous quasars are short-lived active stages of early evolution of galaxies. 

For calculation of bright QSOs number density at different redshifts we have used the 
analytical approach proposed by |32] and improved by [23j. Let us suppose that masses of 
protogalactic clouds (protohalos) which will host brightest QSOs are in the range M, M + AM 
and masses of homogeneous central parts of perturbations (top-hat) in which the SMBH will 
be formed are in the range Mj/^, Mth + AMt/j. The mass of SMBH can be lower than Mth- 

The halo number density with mass in the range Mth-, Mth + ^Mth at some time t is 
determined by the mass function [32j 
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where = a'^{Mth) is the r.m.s. of density perturbations in the top-hat sphere with radius 
R^j^ = (3Mth/4:n'py^^ and 6c(t) is the critical magnitude of the linear density perturbation which 
collapses at t. We will take into account only that haloes with mass Mfh which are central peaks 
of protogalaxy clouds with mass M. They can be extracted from the total number density by 
estimation of the probability that its mass after collapse and virialization of protogalaxy cloud 
will be in the range M, M + AM during short time At |23) 
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where o"^ = a^{M) is the r.m.s. of density perturbations in the sphere which contains the mass 
M of the whole halo. 

The important condition for SMBH formation is spherical symmetry and top-hat density 
profile of central region of protogalaxy cloud in which the appearance of bright quasar is ex- 
pected. All shells in such region collapse practically simultaneously. This should result in 
formation of SMBH. As the main fraction of mass consists of the collisionless cold dark matter 
particles, the pressure of hot baryonic gas can not prevent the collapse. So, from estimated 
number density of haloes with given mass we must extract the essentially non-spherical ones. 
The distribution of the ellipticity e and prolateness p of ellipsoidal clouds with given initial 
amplitude of density perturbation 6 has been obtained by [U [35] : 
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^ ^-2.5(V<x)2(3e2+p2) 



where < e < oo and — e < p < e. We suppose that ellipticity of homogeneous regions of 
protogalactic clouds, the collapse of which leads to SMBH formation, is constrained from above: 
e < eL. Their fraction is as follows 



P(e <e'^) = j de j dpg{e,p\6). 



(2) 



This integral has functional presentation in the formp(e < e^) = /(x), where x = e'^6/a{Mth)- 
It means that constraining of ellipticity of central homogeneous regions e < constrains the 
ellipticity of the whole protogalactic cloud e < e^, where e'^/a{Mth) = em/cr{M). 
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Finally, we present the QSOs number density in the co-moving space as product of the 
number density of haloes with mass in the range Mth-, Mth + AM^/i, the fraction of those of 
them which belong to larger haloes with mass in the range M, M + AM and the fraction of 
such haloes with ellipticity lower than Cm'- 
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We present the dependence of QSOs number density on redshift in the form: 
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All redshift-independent multipliers in are collected by A {[A] =Mpc^^). It is the product of 
relation some constants, spectral dependent values (cr, (Xth), their derivatives with respect 
to mass as well as of the unknown values such as the QSO lifetime {tqso = At) and mass range 
of haloes (AM) which maintain the bright QSOs or X-ray AGNs. Therefore, we assume it is 
a free parameter of model. The critical amplitude of density perturbations 6c{tcoi) depends on 
matter density and value of cosmological constant [21, 20j. For ACDM-model with parameters 
used here it equals 1.673 at 2; = 0. For other redshifts we have approximated the numerical 
results for 6c{z) from j21] by the simple analytic formula 
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which works well in the range of redshifts < 2; < 5. So, the dependence of QSOs number 
density (jlj) on redshift has 3 free parameters: the normalization constant A, mass M and 
maximal ellipticity of the protogalactic clouds e^- They are used for the fitting of ngsoiz) to 
the observational data on QSOs number density. 



2 Results and discussions 

We have used the observational data on bright QSOs number densities at different redshifts 
presented by [36J, which are the compilation of optical (SDSS, COMBO-17) and X-ray (Chan- 
dra, ROSAT) QSOs and AGNs redshift surveys (see FigW)- The ChaMP+CDF+ROSAT data, 
which include the X-ray AGN with Lx > lO^'^'^ erg-s~^ in the energy range 0.3 — 8 keV, are 
shown in the figure by squares. The results of calculations of ngsoiz) for ACDM-model with 
parameters n,, = 0.044, = 0.258, = 0.742, h = 0.719, ag = 0.796 and = 0.963 
[19] are shown there by lines. We calculated the r.m.s. of density perturbation as follows: 
ct(M) = ^(27r)-3/2 k^'^+^'^T'^{k)W^{Rk)dk, where R is the radius of the top-hat sphere con- 
taining the mass M, W{x) = 3(sinx — xcosx)/x^ is the window function and T{k) is the 
transfer function. The latter has been calculated using the publicly available code CAMB 
|24] . The normalization constant of power spectrum A we determined from the equality 
al = >l(27r)-3/2 joo ki^+^')T\k)W\8k)dk. 

The solid line in Fig{T] shows nQso{z) calculated for M = 8 ■ 10^^ h~^MQ and = 0.38. 
The parameter A has been used for fitting of nqsoiz) to observational data at the maximum 
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Figure 1: The dependence of QSOs co-moving number density on redshift: theoretical curves 
versus observational data. Solid line shows the best fitting model with M = 8 ■ 10^^ h'^M^ and 
Cm = 0.38. Dashed and dash-dotted lines present the models with the same protogalactic clouds 
masses M = 4 ■ 10^^ h~^MQ, but different ellipticities, e.^ = 0.2 and Cm = 0.17 respectively. 
Dotted line shows the co-moving space number density of QSOs for the model with M = 
3 ■ 10^1 h-^MQ and = 0.5. 

{z ^ 2.5) and for this line it equals 9.75 • 10^'' Mpc^^. It is the product of lifetime of QSOs, 
spectral dependent values cr{M), a{Mth) and mass range of haloes AM which maintain bright 
QSOs. So, in the ACDM-model the halo of mass M = 8 ■ 10^^ h~^MQ with top-hat central part 
Mth = 10^ H-^Mq may have SMBH with Msmbh ~ 0.1 - O.OOlMi,, which should shine like 
bright QSO during tqso ~ 5.6 • 10^ years (AlogM = Alog Mth = !)• As we can see the solid 
line agrees very well with the observational data ChaMP+CDF+ROSAT |36] covering the wide 
range of redshifts from 2; ^ to 2; ^ 5. 

Without constraint of ellipticity of protogalaxy clouds hosting QSOs the solid line in Fig. 
[1] should be essentially higher than corresponding observational data ai z < 1.2. As it follows 
from ([T]) the mean value of ellipticity of cloud which collapse at given 2 is e = (j{M) / {25c{z)). 
For high z the ellipticity is small as 5c{z) is large. The clouds which collapse at lower z are 
more elliptical and more rarely produce the QSOs. 

The mass M of protogalaxy clouds, which can host bright QSOs, depends on as and 6c{z) 
and should be different in the cosmological models with different sets of parameters. For higher 
(Tg the mass of protogalaxy will be higher. Increasing of the mass leads to the steeper curve 
of QSOs number density ai z > 2.5 and to the displacement of its maximum to lower z. On 
the other hand, the decreasing of the upper limit of ellipticity of protogalaxy clouds leads to 
displacement of that maximum to higher z. The dashed and dash-dotted lines in Fig{T] show 
the QSOs number density for ellipticities e < Cm = 0.2 {A = l ■ 10~^ Mpc~^) and e < = 0.17 
(y4 = 1.1 ■ 10~^ Mpc"'^) correspondingly. The mass of protogalaxy cloud is M = 4 ■ 10^^ h~^MQ 
for both cases. 

The minimal mass of halo in which the bright QSOs may arise can be estimated by com- 
parison of the model curve for quasar number density with corresponding observational data. 
Indeed, the number density of haloes with lower mass is higher and they collapse earlier. So, for 
large z the predicted number density curve will be less steep and its maximum will be shifted 
to z > 2.5. The dotted line shows the redshift dependence of number density of bright quasars 
which arise in the halo with mass M = 3 • 10^^ H^^Mq and initial ellipticity lower than = 0.5 
{A = 5.35- 10~^ Mpc^^). It lays above all observational data (excluding one ROSAT point [29]) 



5 



and that mass can serve as rough estimation of lower hmit for mass of haloes which produce 
the bright QSOs and/or X-ray AGN. 

The dependence of bright QSOs number density on redshift obtained here matches the 
observational data essentially better than ones deduced in the frame of semi-analytical merging 
models of SMBHs formation (see for comparison Fig. 9 from [T8], Fig. 3 from [26], Fig. 7 from 
|27) . Fig. 3 from [5j), which underestimate the number density of luminous AGNs/QSOs at 
high redshift s. 

We must note that number density of quasars with lower luminosity, L2-8keV < 10^^ erg-s~^, 
has maximum at z ~ 1 |2l [HI [151 ISHl HQ]. This can not be explained by simple model for bright 
quasars used here. Some reasons should be taken into account for this case. 1) The assumption 
that lifetime of bright QSOs is small in comparison with cosmological time can be incorrect for 
fainter quasars. It can depend on density profile of protogalaxy cloud, angular momentum etc. 
2) SMBHs can be formed in the protogalaxies of smaller mass essentially later then collapse 
has occured. The time-lag can depend on physical parameters of protogalaxy cloud too. 3) 
Fainter quasars can be the result of reactivation of SMBH in the centra of some galaxies by 
tidal collision of satellite galaxies or their merging. Some galactic nuclei can experience such 
reactivations several times during the galaxy life. 

3 Conclusions 

We have assumed that the bright QSOs are short-lived active phases of early evolution of mas- 
sive galaxies. Their large X-ray and optical luminosities are provided by accretion of baryonic 
matter to SMBHs which have been formed as result of collapse of homogeneous central parts 
of protogalaxy clouds (haloes). Using the expression for co- moving number density of haloes of 
mass Mth present at time t [23] and constraining the ellipticity (e < em) of those of them which 
belong to larger haloes of mass M we have obtained the formula for dependence of QSOs num- 
ber density on redshift. It has 3 free parameters (M, and normalization constant A) which 
allow us to fit the model dependence to observational data. As background cosmological model 
we have used the ACDM-model with best fitting parameters determined by [19j : Q^n = 0.258, 

= 0.742, h = 0.719, as = 0.796 and = 0.963. Obtained dependence for nqsoiz) (HD 
fits in the best way the observational data on bright QSOs and X-ray AGN number densities 
ChaMP+CDF+ROSAT |36] for M = 8 ■ 10^^ H-^Mq and = 0.38. The best fitting value of 
normalization constant is A = 9.75 ■ 10^^ Mpc^'^. The lifetime of such QSOs tqso ~ 6 ■ 10^ 
years. The lower limit of halo mass estimated in such approach is M = 3 ■ 10^^ h~^MQ. 

Therefore, the halo model of galaxy formation together with the assumption that bright 
QSOs are short-lived active phases of early evolution of massive galaxies quantitatively agrees 
well with corresponding observational data. 
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